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ABSTRACT 

We derive a self similar solution for the propagation of an extreme relativistic 
(or Newtonian) radiative spherical blast wave into a surrounding cold medium. 
The solution is obtained under the assumption that the radiation process is fast, 
it takes place only in the vicinity of the shock and that it radiates away a fixed 
fraction of the energy generated by the shock. In the Newtonian regime these 
solutions generalize the Sedov- Taylor adiabatic solution and the pressure-driven 
fully radiative solution. In the extreme relativistic case these solutions generalize 
the Blandford-McKee adiabatic solution. They provide a new fully radiative 
extreme relativistic solution which is different from the Blandford-McKee fully 
radiative relativistic solution. This new solution develops a hot interior which 
causes it to cool faster than previous estimates. We find that the energy of the 
blast wave behaves as a power law of the location of the shock. The power law 
index depends on the fraction of the energy emitted by the shock. We obtain an 
analytic solution for the interior of the blast wave. These new solutions might 
be applicable to the study of GRB afterglow or SNRs. 

Subject headings: Gamma rays:bursts — hydrodynamics — ISM: jets and 
outflows — relativity — shock waves — supernova remnants 



1. Introduction 

Many astrophysical phenomena (SNRs, 7-ray bursts - GRBs, AGN hot spots, etc.) are 
believed to involve radiative shock waves. These shocks accelerate particles which emit the 
observed radiation. In particular, it is widely accepted that the recently discovered GRBs 
afterglow results from an emission by relativistic shocks, created by the interaction between 
an initial ejecta and the interstellar medium. The recent observations of GRB afterglow 
have lead to numerous attempts to model this phenomena. 

If the radiative mechanisms are slow compared to hydrodynamics time scale, the 
blast wave evolution is adiabatic. The propagation of such a blast wave is described be 
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a self similar solution: The Sedov- Taylor solution (|Sedov 1969| , [lay lor 1950|) describes 



the Newtonian regime and the Blandford & McKee (1976) solution describes the extreme 
relativistic regime. 

We call the solution radiative if the radiative mechanisms are fast compared to the 
hydrodynamic time scale. A fully radiative solution is one in which all the energy generated 
by the shock is radiated away. |Ostriker fc McKee ( 1988)1 have shown that if a fully radiative 



blast wave emits its energy only in the vicinity of the shock, it can be described by one of 
two possible self-similar solutions: The pressure driven snow-plow (PDS) or the momentum 
conserving snow-plow (MCS) solution. In the PDS solution a thin shell "snow plows" 
through the external medium, driven by the pressure of its hot, roughly isobaric, interior. In 
the MCS solution the interior has been cooled and a thin shell slows down while conserving 
momentum. [Blandford fc McKee (197B~)| have found an extreme relativistic fully radiative 



solution. This solution describes a thin shell with a cold interior and it can be considered 
as the relativistic generalization of the Newtonian MCS solution. However, momentum is 
not conserved in this solution, as in the relativistic case one has to consider the momentum 
carried by the emitted radiation. 

These solutions are either adiabatic or fully radiative. It is likely that in some cases not 
all the energy would be radiated away, even though the cooling is fast. Such is the situation 
if the shock distributes the internal energy between the electrons and protons and there is 
no coupling between the two afterwards. Since only the electrons cool, only a fraction of 
the internal energy will be radiated. It is likely that at least in the initial phases of a GRB 
afterglow this would be the case. We consider here this "partially radiative" scenario. 

The goal of our paper is to find an analytic solution for the evolution of Newtonian and 
relativistic partially radiative blast waves. This, we believe, will eventually lead to a physical 
description of the evolution of GRB afterglows in their non-adiabatic stage, and of other 
astrophysical phenomena. We find a self similar solution under the assumption that the 
radiative mechanism is fast (compared to the hydrodynamics time scale). We parameterize 
the radiation by a dimensionless (the non dimensionality is essential for a self-similar 
solution) parameter, e, which describes the fraction of the energy produced by the shock 
that is radiated away. We recover the extreme relativistic adiabatic Blandford-McKee 
solution and the Newtonian Sedov- Taylor solution in the corresponding limit when e = 0. 
In the Newtonian case when e — > 1 we recover the PDS solution, but we do not reproduce 
the MCS solution. Similarly in the relativistic limit of e — > 1 the solution approaches a new 
relativistic PDS solution which is different from the radiative Blandford-McKee solution. 



We describe our model in Sec. 0. It is composed of an adiabatic shock followed by a 
narrow radiative region and a wide self-similar adiabatic regime. We proceed in Sec. [O] by 
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calculating the radiative shock conditions for Newtonian blast waves. We describe the self 
similar equations and their solution in Sec. |3.2| . Using the same method we turn in Sec. 
[| to the calculations of relativistic blast waves. Finally, in Sec. |5] we discuss the limiting 
solutions obtained by other authors and compare them with our solution. 



2. The Model 

We consider a spherical radiative blast wave, that appears when a large amount of 
energy is released within a small volume. This results in a strong shock wave that expands 
supersonically into the surrounding medium. We consider the regime where the influence 
of the initial mass is negligible and that the pressure of the surrounding medium is small 
compared to the energy density of the flow. These assumptions are necessary and sufficient 
to obtain a self similar solution, see e.g. |Landau fc Lifshitzt pg. 404. This solution 



corresponds in the adiabatic limit to the well known Newtonian Sedov- Taylor solution or to 
the relativistic Blandford-McKee solution. 

We assume that the shocked, hot matter can radiate a fixed fraction of its internal 
energy generated by the shock. We show that the details of the radiation mechanism are 
not important, as long as the cooling time scale is short compared to the hydrodynamics 
time-scale. We also assume that the radiated energy is not re-absorbed in the system. 
Under these assumptions we can divide the system into three parts. An adiabatic shock, 
followed by a narrow radiative zone and a wide adiabatic flow. 

The assumption of fast cooling allows us to treat the shock and the cooling layer as 
planar and stationary. This means that the velocity of the shock is constant during the time 
that a given fluid element crosses the radiative zone and cools. The physical conditions 
through this "radiative shock" can be found using the equations of energy, momentum and 
particle conservation. For simplicity we assume that the radiation does not alter the shock 
itself, which remains adiabatic, and that a radiative layer follows it. However, the shock 
jump conditions (Rankine-Hugoniot in the Newtonian case and Taub in the relativistic 
case) are derived from the same conservation equations. Thus, even though we divide the 
process into a shock and a radiative layer we effectively use the conservation equations 
between the unshocked matter and the shocked matter after it has been cooled. This is 
valid even if the radiation process changes the shock itself (i.e. in the case where the cooling 
length is comparable to the particles' mean free path within the shock). 

Fig. [I] describes the model. In the shock frame, cold matter enters the shock. The 
shock itself is not affected by radiation and can be described properly by the Newtonian 
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Rankine-Hugoniot or the relativistic Taub jump conditions. The downstream matter is 
hot and it moves sub-sonically relative to the shock front. This matter cools by radiation 
and it is compressed by the surrounding pressure. It stops radiating when it has emitted a 
fixed fraction of its energy. It continues adiabatically afterwards. The assumption of fast 
cooling enables us to treat the shock and the radiation layer as instantaneous and with zero 
width compared to the hydrodynamics scales. This translates to modified radiative jump 
conditions relating the unshocked material and the cooled material after the radiation zone. 
Using these conditions we find a self similar solution for the remaining adiabatic interior. 



3. The Newtonian solution 

3.1. Jump conditions of radiative shocks 

We divide the radiative shock, as defined in Sec. [| to two regions: The "adiabatic" 
shock, and the radiative zone. The conditions for the matter immediately after the shock are 
given by the Rankine-Hugoniot jump conditions. Assuming that the matter is adequately 
described by a polytropic equation of state with an adiabatic index j, the jump conditions 
in the shock frame are (see e.g. Landau fc hifshitz pg. 335) : 



P2=(^\jPu u 2 =(^^U sh) P2 = ^ Pl U^ (1) 

where p\ is the density of the unshocked matter, U s h is the velocity of the shock relative 
to the unshocked matter and p 2 , p 2 and u 2 are the pressure, density and velocity of the 
shocked matter, measured in the shock frame. These jump conditions hold if the shock is 
strong, i.e. the pressure of the unshocked gas is negligible (pi <C p\U^ h ). 

In the radiative zone, the matter that was heated by the shock radiates and cools. In 
the mean time it is also compressed by the surrounding pressure. This compression causes 
the gas to slow down in order to conserve the matter flux. At some point the matter reaches 
some equilibrium and stops radiating. 

To calculate the hydrodynamics of the blast wave interior, we need to know the 
pressure, velocity and density of matter at the end of the radiative layer. Since we assume 
the cooling is fast, the shock velocity is constant during the cooling of a fluid element of 
shocked matter. We look therefore for a steady solution for the conditions of the matter 
during the cooling process. We find the conditions at the end of the radiative layer as a 
function of energy lost through this process. 



In a steady state we can use, in the shock frame, the equations of conservation of 
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matter: 

pu = constant, (2) 

and of momentum: 

(pu)u + p = constant, (3) 

where p, u and p are the fluid's density, velocity and pressure. These two conservation 
equations contain three unknowns, so that we can find the physical quantities as a function 
of one another. We solve equations |^ and |3] for p and p as function of u, 

. , p 2 U 2 

p{u) = , p{u) = p 2 u 2 {u 2 - u) + p 2 , (4) 

u 

where p 2 ,u 2 and p 2 , the density, the velocity and the pressure of the flow immediately after 
the shock, are given by Eq. p]. 

We denote by p, it and p the density, the velocity and the pressure when the fluid has 
stopped radiating, and use the shock conditions ( Eq. [1] ) as the boundary conditions for 
the radiative flow. We use the velocity at the end of the cooling layer to parameterize the 
fraction of energy lost, and find the pressure at that location by 

7-1 2 + £(7-1) Tj2 
p 2 = p 2 {l + —^—5) = -—^ piU sh , (5) 

the velocity relative to the observer, 

U sh -u 2 = (U sh - u 2 ){l + 2^1s) = 2 + 5{l ~ l) U sh , (6) 

2 7 + 1 

and the density is 

P2 7+1 
P2 ~ —5 ~ ( 7 -l)(l-5) Pl ' (7) 

where 5=1 — u 2 /u 2 . Fig. ^| shows the trajectory of a fluid element during the shock and 
the following cooling process. 

During the cooling process energy is radiated away. We write the energy flux as 

J = u(pu 2 /2 + e + p), (8) 

where the energy per unit volume is denoted by e. We find that the fraction of energy flux 
lost in the radiative layer is: 
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Note that we have not used the details of the cooling process to find the pressure and 
density as a function of the velocity. This is sufficient for our calculations, because all we 
need are the physical conditions at a location were the matter has stopped cooling. The 
velocity at this point is given by Eq. [] 

Using Eqs. § and |9] we calculate the overall energy loss rate: 

— = -^R 2 sh J sh e = -2Tr Pl U* h R 2 sh e, (10) 
where J s h = p-JJ^ h /2 is the energy flux of the matter entering the shock in shock frame. 



3.2. The self similar solution 



After the matter has passed the radiation layer, it continues to flow adiabatically. The 
motion is described by matter, momentum and energy conservations equations in spherical 
coordinates: 

dp , 1 9 2 



^_ 1 d , 2 ^ , 
dt r 2 dr 



du du 
dt dr 



0. 



1 dp 

p dr ' 



and 



d_ 

dt 



1 d 



[e + pu 2 /2] + — — [r 2 u(e + p + pu 2 /2)} = 0. 
r z dr 



(11) 
(12) 

(13) 



We look for a self similar solution for these equations. Similarly to the non-radiative 
Sedov- Taylor solution we consider the case where piU 2 h , the momentum flux of the matter 
that enters the shock is much larger than the thermal pressure p\ of the undisturbed 
medium. 

In the Sedov- Taylor solution dimensional arguments lead to a self-similar dimensionless 
variable. In the radiative case the definition of the self similar variable is not straight 
forward, because the energy is not constant in time. We look for a self similar solution to the 
hydrodynamics equations with energy that varies with time as a power-law: E = Eo(t/to) x - 
We define a dimensionless variable 



Pi 



t 2+x {E /t$) 



1/5 



(14) 



1 The cooling profile is determined by the radiation mechanism. It is only the conditions at the end of 
the cooling process that are independent of this mechanism. 



If this variable leads to a self similar solution, r and t will appear in the solution only 
through £, and the position of the shock corresponds to a fixed value of £, denoted £ - Thus: 



R s h{t) = £o 



* 2+a (£q/^ 
pi 



1/5 



dR s h{t) 2 + XR sh 



dt 



t 



(15) 



We substitute R s h and £7^ into Eq. [10] and obtain a cubic equation for A (which has only 
one real solution): 



(16) 



Using the dimensionless similarity variable £, we look for a self similar solution of Eqs. 
of the form: 

p(r,t) =p 2 a(0, (17) 



u(r,t) = (C/ s/l - u 2 ) 
p(r,t) =p 2 ( 



Rsh(t) 



R s h(t) 



(19) 



The coefficients in these equations are chosen in order to match the standard definitions 
in the adiabatic case (see e.g. [Shu 1992|) . Using Eqs. |], §1> we obtain the boundary 
conditions: 

jzr§, «(&) = 1 + 1^5, P (e ) = i + 



a (Co) 



(20) 



To determine the dimensionless position of the shock front, £o ; we require that the 
energy in the blast wave interior equals with the energy defined by the self similar variable: 



E(t) = E \ 



t 



R S h(t) p(r) p(r)u(r) 2 , A 2 , 
I. 7 — 1 2 



(21) 



Substitution of the self-similar functions into this equation yields a non-dimensional 
normalization equation: 



2 + A N9 8tt /■& 



(-£-) 



7^ 



1 Jo 



(p(0 + <M£r)« 



(22) 



Finally we substitute the self-similar functions at(£), t>(£) and into the fluid 
equations 11-13, replacing the partial derivatives with the corresponding total derivatives: 
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We obtain the coupled ordinary differential equations: 
V d log a 



2 d\ogi 

_ d log v 

V — 

d\ogi 

V dlogp 

2"^io g e 



where 



V 



(av(l +j-2v) [(1 + 3A)(1 + j) - 4(2 + X)v) + 
2(A-3)(7 2 -l)p)/(l + 7-2 V ), 

-a(l +j-2v) [5(1 + 7) - 2(2 + X)v] + 
2p(7-l) [(A-3)(1 + 7) + 3(2 + A)t«], 
2(2 + A)(7-l) 7 p + a(l + 7) x 

14 + 2A + 7 + 3A7)w - 5(1 + 7) - 4(2 + A)w 2 
(2(1 - 7 ) 7P + a(l + 7 - 2w) 2 ) (2 + A). 



(24) 

(25) 
(26) 



(27) 

We solve these equations numerically, using the boundary conditions from Eq. and 
obtain the relation between the cooling parameter e and the power-law index A (see Fig. |3|). 
Self similar solutions for several cooling parameters are shown in Fig. 

For e = our solution is adiabatic and we recover the classical Sedov- Taylor solution. 
For e = 1 the density diverges (Eq. 0) and the boundary conditions become singular. Thus, 
to obtain the fully radiative limit we must take the limit of e — > 1. We discuss this limit of 



the Newtonian and the relativistic solutions in section 154. 



For intermediate values of e we see (Fig. ^) that as the cooling parameter becomes 
larger, the matter concentrates in a small shell near the shock and the internal pressure 
decreases. Notice that similarly to the non-radiative Sedov- Taylor solution, the interior 
temperature T oc P/p diverges at the center of the blast wave and it monotonously 
decreases towards the shock. This behavior occurs because the material near the center 
passed through the front earlier, when the shock speed was larger ( infinitely large in the 
formal extrapolation to t — > 0), and it cools slowly, adiabatically. 



4. Extreme relativistic solution 

The model described in Sec. |2| is adequate for relativistic blast waves as well as for 
Newtonian ones. We assume that each shocked particle emits a fraction of the internal 
energy it acquired in the shock. We also assume that the radiation time scale is short 
compared to the hydrodynamic time scale, thus allowing us to treat the blast wave as 
stationary when calculating the radiation process. 
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We start by calculating the conditions at the end of the "radiative shock" as a function 
of the fraction of energy emitted. Then we use these results as boundary conditions for the 
blast wave, assuming that the flow outside of the radiation layer is adiabatic. Following 
B landlord fc McKee (1971)) we derive a self similar solution for extreme relativistic blast 



wave, accurate to the leading order in T 2 . However, within the derivation the second order 
terms in the equations are also important. 



4.1. Jump conditions of radiative shocks 

We begin by calculating the jump conditions for a radiative extreme relativistic strong 
shock (r ^> 1) using the relativistic equation of state 7 = 4/3. The jump conditions are 
derived from the continuity of the energy, momentum and particle flux densities in the 
shock frame, (see e.g. [Landau fc Lifshitz| pg. 511 ). Assuming that the unshocked matter 



is cold, pi <C P1T 2 , we have 

e 2 = 2T 2 Pl , (28) 
p 2 = 2V2 Pl T, (29) 

ll = \? 2 - (30) 

For generality we allow for any polytropic index, and relax our assumption of extreme 
relativistic flow when calculating the cooling profile. Later on, while tying together the 
shock and the cooling layer, we return to an extreme relativistic motion. 

In a relativistic flow the thermal energy can be comparable to the rest mass energy, 
and we have to take into account the momentum of the radiation. We cannot calculate 
the profile from matter and momentum flux conservation alone as in the Newtonian case, 
and the independence of the cooling profile from the nature of the radiation mechanism is 
not obvious. We assume that the radiation is emitted isotropically in the fluid rest frame 
with a local cooling function C(e,p) (like in the Newtonian case this function does not 
appear in the final result), and calculate the cooling profile in the shock frame using matter 
conservation: 

= 0, (31) 

momentum conservation: 

^[7 S2 /5 S2 (P + ie) + (7 - l)e] = -2^W (32) 
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and energy conservation: 



d_ r 

dx 



Y' 2 P S (P + ie)} = -^^l s \l + (3 s2 ). (33) 



To solve these equations we rewrite Eq. pll as 



P = (34) 



The cooling function C can be eliminated from Eqs. [32| and [33], which after a little algebra 
and the usage of Eq. p4| results in: 



d(3 s ~ 7-l-/3 s2 ' 
Integrating this differential equation we obtain 



(35) 



e 2 + , ^ 7 ^ 2 log mP s )/m)\ , (36) 
2^/(2-7X7-1) 



where 



a/ax = (1 - WtHW + Iv 7 ^! + V2^7/7) 
lPJ (l + /g /V7^T)(l-/5 v / 7^T+v / 2^7/7)' 1 J 



Fig. [5] depicts the flow of matter through the "radiative shock" . In the shock frame 
the dense shocked matter flows out of the shock with velocity c/3, slows and becomes even 
denser during the cooling process. In the extreme relativistic limit the pressure does not 
change during this process. 

Combining the radiative flow solution with the strong shock jump conditions we obtain: 

e 2 = 21*01, (38) 

p 2 = 2V2pir^±I, (39) 
7 2 2 = ^r 2 (l + e). (40) 



Using Eqs. S8L p9 and 40 we find that e is equal to the fraction of energy that each 



particle has lost in the unshocked fluid frame during the cooling process: 

e _ 47 2 e 2 /3p 2 ~ 47 2 e 2 /3p 2 
47 2 e 2 /3p 2 



- 11 - 



This definition of e coincides with the definition of Sari (1997) 



The energy loss rate of the whole blast wave in the unshocked matter rest frame is the 
difference between the work done on the cooling layer by the internal pressure p 2 , and the 
increase in internal and kinetic energy of shocked matter. For spherical radiative blast wave 
we obtain: 



dt sh 



P2P2 - {(e 2 +P2W2 -P2] [m -fa}] = -87ri4r 2 e/3. (42) 



This rate equals the rate in which energy is supplied to the unshocked matter (47i\R 2 p 2 ), 
multiplied by e. [Blandford fc McKee [T976J have obtained an equation for dE/dt in the 
limit e — > 1 (Eq. 84 there) which differs from our result due to a missing factor of 4/3 in 
their Lorentz transformation of the energy density. 

It is important to note that there are many possible definitions of radiative efficiencies 
which are frame dependent. For instance, the fraction of internal energy lost in the shock 
frame is: 

e 2 (4 7 f - l)/3p 272 s - S 2 (4 7 ~f 2 - l)/3p 2 7 2 s e(37 + e(16 + 3e)) 



e 2 (4 7 | 2 - l)/3p2 7 | " 7(l + e)(3 + e) 

This equals e for e = and e = 1, but it has a different form for intermediate values. 



(43) 



A physical example of a partially radiative shock wave leads to another definition of 
radiative efficiency Consider a system in which the shock distribute the energy between 
electrons and protons according to: 

ef c = e e e 2 , e p 2 rot = (1 ~ e e )e 2 , (44) 

and there is no coupling between electrons and protons afterwards. Charge neutrality 
requires that the density of electrons and protons remains equal. The protons do not 
radiate. Therefore their flow is adiabatic. When the electrons have cooled down completely, 
the protons' energy satisfy: 

= e Vrot p -^ (45) 

At this stage the cold electrons have no internal energy, e|' ec = and therefore e 2 = e 2 ro *. 
Using Eq. [39| we obtain: 

«■ = 1 " fe)' (46) 



2 Energy density transforms like e' = (e+p)~f 2 —p, which becomes e' = (4/3)e7 2 in the extreme relativistic 
limit. In this calculation Blandford fc McKee (1976) have used e' = ej 2 . 
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Again, e e equals e for e = and for e = 1 (see Fig. |9|), but it has a different form. It is simple 
to understand this by recalling the behavior in the cooling layer. As a fluid element cools 
it slows in the shock frame and accelerates in the observer frame. The electrons accelerate 
the protons while cooling in order to maintain charge neutrality. This acceleration transfers 
energy from the electrons to the protons. Therefore, not all the electrons internal energy is 
radiated away. 



4.2. Self similar solution for homogeneous medium 

Inside the blast wave the flow is adiabatic and the equations of motion are obtained by 
setting the divergence of the energy momentum tensor to zero (see e.g. |Landau fc Lifshitz| 
pg. 506 ). For spherically symmetric flow, the pressure, energy and velocity satisfy 

W + ^d~r r ^ = ° (47) 
d 1 d 

- 7 2 (e + p + p)P + -^^{e + p + p){3 2 = (48) 
d 1 d 

- 7 2 (e + p + (3 2 p) + ^r 2 7 2 (e + p + p)/3 = (49) 

where (r, t) are in the unshocked matter frame and r = is the center of the blast wave. 

As in the Newtonian case we look for a solution in which the energy decreases as a 
power-law in time: 

T 2 oc r m , m > 3. (50) 
Following [Blandford fc McKee (197571, we define a similarity variable 



X = [l + 2(m + l)r 2 (t)](l-r/t). (51) 

We substitute the self similar variables and take the extreme relativistic limit f\ of Eqs. 
47L [4~8| and [49]. We have to expand Eqs. [48| and E9| to the second order in T -2 as these 



equations are identical in the leading order in T 



-2 



3 While solving Eqs. |^] and in the extreme relativistic limit we assume p <§; e. We check the validity 
of this limit by requiring that the flow continues to be extreme relativistic after the cooling layer, e^jpi 3> 1. 
This results in a requirement r>l/(l-6). It means that as the blast wave looses more energy, our solution 
breaks down, even though f>l. 
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We write the pressure, velocity, and density of the shocked fluid as 

p(r,t) 7 (r,t) = %p 2 h( X ) = 2 Pl T{t) 2 \±^h{ X ). 

1 — e 



(52) 

(53) 
(54) 



These expressions are valid for the interior of the blast wave, x ^ 1- The boundary 
conditions at X — 1 (Eqs. [58^ ^ and [40]) are satisfied by setting 

/(l) = 0(1) = /!(!) = 1. (55) 



Substituting Eqs. [52], |53] and |54| in Eqs. [47], [48] and [49], and replacing the derivatives 



using 



t 



dT(t) mT{t) 
dt ~ ~^2~T' 

= (m-L 1 V9F 2 

dt 



(m + l)(2r 2 - X ) + l, 



and 



^ = -[l + 2(m + l)r 2 ] 



all to the second order in (r 2 ), we obtain the following self similar equations: 
A ding 



9 dx 
A din/ 

9 dx 
Adlnh 

9 dx 



where 



A 



= (7m - 4) - (m + 2)(1 + e)gx 

= 8(m- 1) - (m-4)(l + e)#x 

= (2(9m-8) -2(5m-6)(l + e)#x + 
(m-2)(l + e )Vx 2 )/(2-(l + e)^x), 

(m + 1)(4 - 8(1 + e)^x + (1 + e) W) 
(1 + e) 



(56) 
(57) 

(58) 

(59) 
(60) 
(61) 



(62) 



To solve these equations we need a relation between the amount of cooling, presented 
in the equations by the parameter e and the evolution parameter m. The energy stored in 
the blast wave is given by: 



E(t) 



— r e 4 7 r,t - l)dr = —, — v 

3 3(m + 1) 



r¥(i + e) 



fgdx, 



(63) 
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to the leading order in T . The time derivative of the energy satisfies 

dE E 

- = (3- m ) T . (64) 

Combining Eqs. ^2], ^ and we obtain a normalization equation which combines 
the hydrodynamic profile through / and g, the evolution parameter m, and the energy loss 
parameter: 

i- 7 ^ Tfgdx-T^-. (65) 

m + 1 Ji 1 + e 



Surprisingly, there is an analytic solution to Eqs. |6(], |61] and |65] for arbitrary values 

of e: 

e 2 + 14e + 9 

m = (66) 

9 = X~\ f = X' a \ h = X ~ a \ (67) 

5 e 2 + 14e-21 

Oil — 1 H , CKo = (68) 

12 + e' e 2 + lle-12 v ; 

and 

E(t) = ^Pl^±YH\ (69) 

W 3 17 + e V ; 

The functions /, g, h are displayed in Fig. ^. In the limit of an adiabatic blast wave, e — > 0, 

we recover the Blandford-McKee solutions. The matter is concentrated, in this solution, 

in a narrow shell of width R/T 2 . As e increases, the matter tends to form even narrower, 

denser shell. Similarly to the Newtonian solution, the internal energy does not decrease 

as fast as the matter density in the interior of the blast wave. In fact, the width of the 

pressure profile is approximately independent of e. 

Our solutions are expressed in the un-shocked matter frame, which is the same as the 
observer frame. However, due to the relativistic motion of the emitting matter toward the 
observer, the time difference between the arrival of two photons to an observer at infinity 
is not equal to the delay between their emission. A photon released at time t from the 
radiation layer will reach the observer at a time 

i * = i - fl < t > = 2(^TT)F 0Ct '" +1 - < 70 > 

All the observables (luminosity, frequency, etc.) of blast waves should, of course, be 
given in the observer time. For instance, the bolometric luminosity is the derivative of the 
blast wave energy to the observer time. Using Eq. |64] we find 

dE -^if-i 
Loc-— ~t £ +1 . (71) 

^obs 
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An interesting case is a sudden release of all the blast wave energy. As the energy is 
released, the shell Lorentz factor also drops to unity. Due to the dependence of the observed 
time on T, this event, no matter how short it is, will be spread over a long time. Even if 
m — > oo, the observed luminosity will not drop faster than L oc t~ 2 s . (In our solutions this 
limit can not be reached because m < 12 for all e). In Fig. [7] we show the relation between 
the energy drop rate and e. 

We check our solution for the blast wave interior using a spherical one dimensional 
numerical simulation. For simplicity, to avoid detailed modeling of the cooling process, 
we use our analytic solution for the "radiative shock" . At each time step we calculate the 
location we expect the shock to be, assuming that the evolution follows our self similar 
solution. We set the values of density, energy and velocity at that location using the 
modified shock conditions of Eqs. |3S|-f|U| As the initial conditions we use our self similar 
solution, and check if the calculated profiles follow it. We continue until the solution 
becomes non extreme relativistic. We find a good agreement between the numerical profiles 
and the self similar ones, see Fig. |8|. 

Blandford & McKee (1976) have also found a self similar solution for blast waves 
with injection of energy. This solution contains an additional internal shock wave. It is 
interesting to check whether it is possible to find another radiative solution by incorporating 
such a shock wave into our self similar solution. To check this, we write the velocity of 
a sphere with constant x by inverting Eq. |5ll to r — r(x,t) and applying a partial time 
derivative. We obtain: 

it = r 2 /x- (72) 

The velocity difference between the fluid and this sphere is: 

r<«// = — , , (73) 

2J2 1 + 6 



which is constant over the whole profile. We find that for all possible values of e, the 



Lorentz factor of the velocity difference Tdiff < y 3/2 (the local speed of sound), i.e. the 
flow is subsonic, while a shock requires supersonic flow. Therefore an additional shock wave 
cannot take place within this self similar solution. 



4.3. External medium with a density gradient 



These solutions, described in Sec. [4.2| , can be generalized to the case where the external 
medium density has a power law profile. If the density gradient is p% oc r~ k , we have 

pi oc T 2k/m , E oc r 2+2k/m R 3 oc iT (m - (3_fc)) . (74) 
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The self similar are: 

A ding 

9 dx 
A din/ 

9 dx 
Adhih 

9 dx 



(7m + 3k-4) - (m + 2)(l+e)gx (75) 

8{m-l)+Ak-{m + k-A)(l + e)gx (76) 

(2(9m + 5k - 8) - 2(5m + Ak - 6)(1 + e)#x + (77) 

(m + fc-2)(l + e) W)/(2 - (1 + e)^), (78) 



where A is defined by Eq. [62|. Equations |f5H78| have an analytic solution: 

9 = X~\ f = X~ a \ h = X ~ a2 (79) 

17 + e-4A; 21 - e(14 + e - Ak) - 6k . . 

ai= 12 + e-3A;' " 2 = (1 - e )(12 + e - 3Jfc) (80) 

for 

m= (l + t )» + 3(l + 0(4-t)-4 > _ L 
3 — e 

The energy contained in the blast wave is therefore given by: 

E(t) = g^i 3 ~ e - r¥, (82) 

w 3 17 + 6-4A; v ; 

where pi should be understood as the density of the external medium at the position of the 
shock. 

As in the solution with homogeneous medium, it is interesting to check if another 
solution can be obtained by incorporating another shock into the self similar solution. 
However, Eq. [7^ is valid also for the this solution, and no additional shock can be fitted 
into this self-similar solution. 



5. Comparison to the limiting radiative cases 

5.1. The fully radiative solution 

The fully radiative solution deserves a special attention due to the singularity in the 
boundary conditions both in the Newtonian and in the extreme relativistic limits. To clarify 
the situation we re-derive following |Ostriker fc McKee (1988") the Newtonian fully radiative 



solutions. Then we discuss the validity of these solutions and compare them to our self 
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similar solution. Next we obtain a new fully radiative solution for the relativistic case and 
compare it to the Blandford-McKee fully radiative solution and to our self similar solutions. 

In the Newtonian case we assume that all the matter is concentrated in a narrow shell, 
and look for a self similar solution. We define a dimensionless kinetic energy, and use Eq. 



16 with e = 1 to find 

a - WW - -T<-—x h (83) 

where V = AirB? sh /3 is the volume swept by the blast wave. The pressure must be constant 
in the empty interior (otherwise it will lead to infinite acceleration). We find 

P = ) A - \piVUl) = Pl Ul h {- - i)(7 - 1). (84) 

V to Z a Z 



The equation of motion for a narrow shell is: 

d( Pl VU sh ) 



4nR 2 sh P. (85) 



dt 

Substitution of Eqs. |3| and |84] in Eq. ^5| yields a quadratic equation with two solutions: 

x J "I MCS <*fi\ 
PDS (86) 

The MCS solution incorporate an empty interior, with no pressure. In the PDS solution 

during the blast wave expansion corresponding to an adiabatic expansion in which new 
matter does not enter the blast wave interior. This behavior is consistent with our 
assumption that new matter accretes only on the expanding shell. 

The existence of two possible solutions with the same boundary conditions is unusual. 
Therefore, one should recognize the physical conditions which leads to each of the self 
similar solutions. |Cioffi et al. ( 19887] have found numerically that a supernova remnant 
(SNR) evolves initially according to the adiabatic Sedov- Taylor solution, becomes PDS 
when the energy loss near the shock becomes important, and finally becomes MCS when 



other processes cool the blast wave interior. Gaffct (1983) has investigated blast waves with 



power-law cooling functions, C oc pT~ c . He has found that c > —2/3 leads to the MCS 
solution, c — > —2 leads to the PDS solution and other values of c do not lead to self similar 
solutions. Since lower values of c correspond to less effective cooling in lower temperature 
(i.e., to a more adiabatic interior), no interior cooling results, asymptotically, in the PDS 
solution. 
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Approaching the limit e — > 1, without taking into account interior cooling, results in 
the PDS solution. During the PDS expansion no matter enters the interior which expands 
adiabatically. Incorporating even an infinitesimal interior cooling will cool the interior down 
eventually, and the blast wave will asymptotically approach the MCS solution. Therefore, 
the limit of a fully radiative blast wave with adiabatic interior results in the MCS or PDS 
solution depending on the way we reach the limit of an adiabatic interior.. Our self similar 
solution reaches the PDS solution in the e — > 1 limit. This is due to the adiabatic interior 
assumption, which, following the previous arguments, forbids a MCS solution. 

Blandford & McKee (1976)' have calculated the dynamics of a fully radiative relativistic 
blast wave with an impulsive supply of energy, treating the swept-up matter as lying in 
a thin, cold shell adjacent to the shock and assuming a cold interior (similarly to the 
Newtonian MCS solution). They have found that T cx -R -3 , which translates to an evolution 
of energy as seen by the observer: 

E oc tJJ 7 . (88) 

We proceed by deriving a PDS like relativistic solution. In the extreme relativistic limit 
we assume that e ^> p in order to obtain self-similar solution. The motion of a matter shell 
is not self-similar, and we cannot follow the Newtonian derivation step by step. However, 
we know from Eq. [73] that in the e — 1 limit the matter moves with same speed as shell 
of constant Xi i- e - new matter does not enter the blast wave interior, as in the Newtonian 
solution. We obtain the internal pressure, P, looking at the evolution of a shell with a 
constant x an d a constant width d\- Using Eq. [TTJwe find to the leading order in T~ 2 : 

dR oc d X T- 2 - 2/m , R oc r~ 2/m , (89) 

where dR and R are the shells' width and its radius, in the unshocked matter frame. We 
use Eq. [52] to find that P oc T 2 , and obtain: 

PV* oc r 2 ~ (6/m+1 ^, (90) 

including a Lorentz transformation of the shell width to the matter frame. Requiring that 
PV 1 = const, we obtain m = 12. Using Eq. [fl] we finally find: 

E^tJJ 13 , (91) 

As in the Newtonian case, our self similar relativistic solution approaches the PDS 
solution in the fully radiative limit. An analytic treatment of the validity of the MCS and 
PDS solution does not exist in the relativistic case. However, in view of the Newtonian 
solution we expect that interior cooling is needed in order to obtain the MCS solution. 
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5.2. Solutions for e C 1 

Using the conditions at the shock, and our radiation model with parameter e, we have 
found the energy loss rate for the Newtonian (Eq. pH| ) and for the extreme relativistic (Eq. 
42) limits. Assuming that a self similar solution exists, we have found scaling laws for the 



conditions at the shock, (Eq. [y| or Eq. |50]). Substituting the energy (calculated from the 
self similar profiles, Eqs ^1] or |63] ) and combining the two, we have found equations which 
relates the energy decrease rate, the radiative parameter and the hydrodynamic profiles 
(Eqs. [16] or Eq. 



Self similarity requires that the energy is a power-law of time. The non-dimensional 
profiles are needed to determine the power-law index. In Sec. EO and Sec. (O of this paper, 



we have solved the self similar equations to obtain this index. However, as suggested by the 
relativistic calculations of pari (1997)1 , if e 1 a simpler derivation exists. We assume that 



in this case the radiation does not alter the hydrodynamic profiles and use the adiabatic 
Sedov- Taylor or Blandford-McKee profiles to obtain the power-law index. 

In the Newtonian case the hydrodynamic profiles enter Eq. [16] through the self similar 
location of the shock front. Assuming that the profiles do not depend on radiation, the 



shocks' location is constant, i.e. £o = £o(e = 0). Substituting this into Eq. [16] we find the 
energy decrease rate to the lowest order 

For 7 = 5/3 we obtain A = — 0.81e + 0(e 2 ). The limiting curves appear in Fig. |[ 

In the relativistic case the profiles enters Eq. ^5] through f£° j gd\. We use the 
adiabatic Blandford-McKee solution to obtain J£° fgd\ = 12/17, and find to the lowest 
order in e 



, x , , e + 0(e 2 ) = — e + 0(e 2 ). (93) 
dlog(t obs ) J^fgdx 1 ; 12 1 ; 1 ; 

This limit is shown in Fig. [7[ This result differs from [Sari (199T)| due to a missing factor of 



4/3 in his Lorentz transformation of the energy density (see footnote [|). 

We conclude by noting that in the Newtonian and in the extreme relativistic limits, 
the change in the blast wave evolution causes the energy to decrease slower than what was 
calculated earlier using adiabatic estimates. 
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6. Conclusions 

We have solved the hydrodynamical evolution of blast waves in which each shocked 
particle emits a fixed fraction of the energy it gains in the shock. We have divided the blast 
wave into three regions: the adiabatic shock, the radiative layer, and the adiabatic interior. 
The adiabatic shock and the radiation layer were combined to form a "radiative shock" , 
which set the boundary conditions for the adiabatic interior. 

For fast cooling cases we have obtained a solution for a planar radiative layer with 
arbitrary shock velocities, independent of the cooling process. We have found that in 
the shock frame the fluid cools, slows down and its pressure increases during the cooling 
process. 

We have obtained self similar solutions for adiabatic interior for the Newtonian and 
the extreme relativistic cases. We find that radiation can change the hydrodynamics 
considerably. Because of this change in the blast wave evolution the luminosity decays 
slower (in time) than what was estimated earlier assuming that the self similar profiles are 
independent of the radiated energy. 

We have also found that as a blast wave becomes more radiative, its matter concentrates 
near the shock and forms a dense shell. However, the internal pressure does not drop to zero. 
In the fully radiative limit of Newtonian blast waves we have reached the pressure-driven 
solution. We have obtained a new extreme relativistic solution for fully radiative blast 
wave, which resembles the Newtonian PDS solution. Our self similar solutions reach this 
modified solution in the fully radiative limit. This solution does not correspond to the 
Blandford-McKee radiative solution. 

The pressure must be continuous within the adiabatic interior and it cannot develop 
self a similar shock or a rarefraction wave (apart from the main strong shock with the ISM). 
Therefore, even in the Newtonian PDS limit with isobaric interior, the solution contains 
a self similar transition layer between the "radiative shock" conditions and the internal 
pressure. 

The recently discovered X-ray, optical and radio emission following a GRB, so called 
"Afterglow" is widely believed to be the result of the deceleration of a relativistic material 
that collides with the surrounding matter. According to the common model, the shock wave 
produced by the collision accelerates electrons to relativistic velocities. These electrons, 
that carry a fixed fraction of the internal energy produced by the shock emit synchrotron 
radiation which is the observed afterglow. For reasonable parameters, the electron cooling 
is fast (at least in the early stages of the afterglow evolution), so that the electrons loose 
most of their energy. Since the density behind the shock is low, the coupling between 
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the electrons and the protons is negligible. Therefore, the protons energy is not radiated 
away. This leads to a partially radiative rather than a fully radiative shock, followed by 
an adiabatic flow. This is exactly the scenario that leads to the partially radiative blast 
waves derived in this paper. Our new self similar partially radiative blast wave can serve 
as the basic hydrodynamic solution from which spectra and light curves of the afterglow 
can be calculated. In particular we derive a new T(R) relation (Eqs. |66|, |5(] ) between the 
Lorentz factor and the radius. This is the critical relation that determines most afterglow 
observations. 

This research was supported by a US-Israel BSF grant 95-238 and by a NASA grant 
NAG5-3516. Re'em Sari thanks the Clore Foundations for support. 
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Fig. 1. — Density as a function of distance from the center of the blast wave. The schematic 
drawing depicts the model used to obtain a self similar solution. It includes an adiabatic 
shock, a small radiative regime and an adiabatic self similar regime. 
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Fig. 2. — Pressure as a function of density for the shock and the cooling profile. The dotted 
arrow illustrates the adiabatic shock. The fluid then follows the appropriate line ( Newtonian 
- thin; Extreme relativistic - thick) until it reaches the cooling parameter. Labels A and B 
describes the density and pressure of the shocked matter immediately after the shock (p2, 
P2)- Note that in the relativistic case the pressure is constant within the cooling layer. Label 
C describes the maximal pressure (piU] h ), reached in the fully radiative Newtonian case at 
the end of the cooling layer. 
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Fig. 3. — The energy loss rate —d\og(E)/d\og(t) as a function of e for 7 = 4/3 (lower 
curve) and 7 = 5/3 (upper curve). The dashed-dotted lines are the linear approximation 
for e — > 0, assuming that the self similar profile does not depend on e. The circles are the 
corresponding PDS solutions. 
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Fig. 4. — The normalized density (solid), velocity (dashed-dotted) and pressure (dashed) 
profiles of the self similar Newtonian solution, for e = 0, 0.3, 0.7, 0.9. The arrows show the 
direction of increasing e. 
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Fig. 5. — The velocity in the shock frame as a function of the density for the extreme 
relativistic shock (dotted arrow) and for the cooling profile (solid line) 
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Fig. 6. — The normalized density (solid), the Lorenz factor (dashed-dotted) and the pressure 
(dashed) profiles of the self similar extreme relativistic solution, for e = 0, 0.5, 0.8, 0.95. The 
arrows shows the direction of increasing e. 
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Fig. 7. — The observed energy loss rate, —dlog(E)/d\og(t t, s ), as a function of e, extreme 
relativistic case (solid line). The dashed line is the energy loss rate for the fully radiative, 
MCS like, Blandford-McKee solution. The dashed-dotted lines are the linear approximation 
for e — > 0, assuming that the self similar profile does not depend on e 
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Fig. 8. — Simulated (x-marks) and self similar (solid line) solutions for the interior of an 
adiabatic (e = 0) and a partially radiative (e = 0.8) blast waves. The initial conditions are 
r = y/2- 500, R—l. Note the larger Lorentz factors and densities in the radiative solution. 
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Fig. 9. — Radiated energy as a fraction of the work done by the shock on the surrounding 
medium, vs. the fraction of internal energy the shock distributes to electrons. The calculation 
assume that electrons and protons are not coupled outside of the shock, and that only the 
electrons radiate. 



